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functions (PDFs) of the proton during 2010. We discuss a recently released CT10 family of 
PDFs with a fixed and variable QCD coupling strength; implementation of combined HERA 
and Tevatron lepton asymmetry data sets; theoretical issues associated with the analysis of 
W charge asymmetry in PDF fits; PDFs for leading-order shower programs; and constraints 
on new color-octet fermions from the hadronic data. 
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CTEQ PDFs in 2010 



Parton distribution functions (PDFs) are essential nonperturbative functions of quantum chromo- 
dynamics (QCD). They describe the internal structure of the proton in high-energy scattering at 
the Fermilab Tevatron collider, CERN Large Hadron Collider, and in other experiments. Mod- 
ern PDFs continuously evolve to include emerging theoretical developments and latest data from 
hadronic experiments, and to provide reliable estimates of uncertainties associated with various 
experimental and theoretical inputs. In this paper, we review the recent progress in the determi- 
nation of the PDFs by CTEQ collaboration [1, 2, 3, 4], which is one of three groups involved in 
the global analysis of hadronic data, besides the MSTW [5] and NNPDF [6] groups. 

Implementation of new data sets 

Since the release of the previous general-purpose CTEQ6.6 PDF set [7] in 2008, new data sets 
have been published in every category of processes included in the global QCD analysis: deep 
inelastic scattering (DIS), vector boson production, and inclusive jet production. These data sets 
include a combination of DIS cross sections by the HI and ZEUS collaborations in HERA-1 [8], 
as well as measurements of W lepton asymmetry [9, 10, 11], Z rapidity distributions [12, 13], and 
single-inclusive jet cross sections [14, 15] by CDF and D0 collaborations at the Tevatron. All these 
new data are included in our latest global analysis, designated as CT10 [3]. 

The new analysis produced two families of general-purpose PDF sets, denoted as CT10 and 
CT10W, which differ in their treatment of the Tevatron W lepton asymmetry data sets affecting 
the ratio of d and u quark PDFs at x > 0.1, as discussed below, but are very similar in all other 
aspects. In addition, we examined the dependence of the PDFs on the QCD coupling a s (Mz) and 
provided special CT10AS PDF sets with a varied a s (Mz) in the range 0.113-0.123 to evaluate the 
combined PDF-a s uncertainty in practical applications. The CT10 PDFs are obtained at next-to- 
leading order in a Sl using the general- mass treatment of charm and bottom quark contributions 
to hadronic observables. To support calculations for heavy-quark production in the fixed-flavor- 
number factorization scheme, we also provide additional PDF sets CT10(W).3F and CT10(W).4F, 
obtained from the best-fit CT10.00 and CT10W PDF sets by QCD evolution with three and four 
active quark flavors. All the PDF sets discussed in this paper (CT10, CT10W, CT10AS, CT10XF, 
and CT09MC) are available as a part of the LHAPDF library [16] and from our website [17]. 

Constraints from combined DIS data by HERA-1 

The CT10/CT10W fits include a combined set of HERA-1 cross sections on neutral-current and 
charged-current DIS [8], which replaces 11 separate HERA-1 data sets used in CTEQ6.6 and 
earlier fits. In the combined set, systematic factors that are in common to both experiments were 
presented as a table of 114 correlated systematic errors, whose effect is shared by each data point 
in all scattering channels. As a result of the cross calibration of detection parameters between 
the HI and ZEUS experiments, the combined data set has a reduced total systematic uncertainty. 
Consequently, the PDF uncertainties at x < 10~ 3 , in the region where the HERA data provide 
tightest constraints on the gluon and heavy-quark PDFs, are also reduced. 

The impact of the combined HERA-1 set on the PDFs is illustrated by Fig. 1, showing relative 
differences between the CT10 PDF set, fitted to the combined HERA-1 data, and a counterpart 
fit, fitted to the separate HERA-1 data sets. In the left subfigure, comparing the best-fit PDFs in 
the two fits, one observes reduction in the gluon and charm PDFs at x < 0.05, accompanied by 
a few-percent increase in the u and d quark PDFs in the same x region. The strange quark PDF 
shows a larger suppression (up to 25% at x = 10~ 5 ), which, however, is small compared to the large 



Figure 1: Left: ratios of CT10 central PDFs fitted to the combined HERA-1 data set and to the separate 
HERA-f data sets, at scale jj, = 2 GeV. Right: bands of the PDF uncertainty (relative to the the central 
PDF set) for the gluon PDF in (red) CT10 with the combined HERA-1 data, (blue) CT10 with the separate 
HERA data. g{x), Q = 2 GeV. 
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PDF uncertainty associated with this flavor. Fig. 1 (right) shows the asymmetric fractional PDF 
uncertainty, computed as in [18], and normalized to the best-fit gluon PDF of each fit. The impact 
of the HERA-1 data on the uncertainties of the gluon and charm PDFs is visible in the small- a; 
region, starting from x — 10~ 3 and going down to x = 10~ 5 , where the error bands contract upon 
the combination of the HERA data sets. In the large x region, the error bands for the combined 
and separate HERA data sets are almost coincident. 

Agreement with the HERA data at small and large x 

The overall agreement of the CT10 fit with the combined HERA-1 data is slightly worse than with 
the separate HERA-1 data sets, as a consequence of some increase in x 2 /d.o.f. for the neutral- 
current DIS data at x < 0.001 and x > 0.1. While the origin of this increase is uncertain, the 
pattern of point-to-point contributions to \ 2 from the data is consistent with random fluctuations 
that turn out to be larger than normally expected. No systematic discrepancies between the 
HERA-1 DIS data and theoretical cross sections are observed, suggesting that the NLO QCD 
theory based on CT10 PDFs is generally consistent with the HERA experiments in the region 
Q>2 GeV included in the CT10 fit. 

In looking for potential systematic deviations of this kind, we examined the agreement with the 
data as a function of either x and Q, or the "geometric scaling" variable A gs — Q 2 x 3 proposed by 
NNPDF authors in Refs. [19, 20]. At A gs — > 0, DGLAP factorization that is required to introduce 
the PDFs can be invalidated by higher-twist terms or saturation; the question is whether such 
effects may bias the determination of the PDFs at A gs > 0.1, in the kinematical region commonly 
included in the global fits. 

Indeed, the NNPDF study finds that the PDFs fitted to the HERA data above some cutoff 
value (A gs > A cut ), disagree at the 2a level with the HERA data in the "causally connected" 
region below the cutoff, 0.5 < A gs < A cut . The NNPDF analysis is realized in the zero-mass 
approximation and includes DIS data in the less safe region V2 GeV < Q < 2 GeV. We repeated 
the A cut fits proposed by NNPDF as closely as possible, in the general-mass factorization scheme, 
and in the region of Q > 2 GeV where our data are customarily selected to suppress higher-order 



Figure 2: Comparison of CDF Run-2 lepton asymmetry data [9] with LO, NLO, and resummed 
predictions from ResBos [35]. 
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and higher-twist terms. While the outcomes of our A cut fits bear some similarity to those by 
NNPDF, the discrepancies between our best-fit NLO predictions and the data below A cut are less 
significant than those quoted by NNPDF and are characterized by a large PDF uncertainty. Thus, 
our fits do not corroborate the existence of stable deviations of the NLO DGLAP factorization 
from the data, if the lower Q bound is chosen to be above 2 GeV. See further discussion in the 
appendix of Ref. [3]. 

Tevatron Run-2 W lepton asymmetry data 

The puzzle of Run-2 W asymmetry. Recently, the Fermilab DO Collaboration [10, 11] 
published measurements of W charge asymmetry At{yt) in electron (£ = e) and muon (£ = fi) decay 
channels, presented as a function of the rapidity of the charged decay lepton. NLO predictions 
based on CTEQ6.1 and CTEQ6.6 sets disagree with these data at surprisingly large x 2 /Npb of 
about 5. The values of x 2 /Npt can be even higher (as high as 20) for some other recent (N)NLO 
PDF sets [11, 22]. Such level of disagreement may appear surprising, given that the Tevatron W 
asymmetry probes the ratio of d and u quark PDFs [21] in the region x > 0.1, where they are 
known quite well from the other experiments. 

Sensitivity to the d/u slope. The discrepancy involving At can be understood in part by 
noticing that the At measurement is very sensitive to the average x derivative (slope) of the ratio 
of the up and down quark PDFs, d(x,Mw)/'u(x,Mw), computed between the typical x values 
%i,2 = Mwe ±Vlv I ' y/s accessible at a given boson rapidity yw [21, 23]. Small variations of the d/u 
slopes in distinct PDF sets can change the behavior of At by large amounts [24] . 

Impact of soft gluon resummation. Another factor at play are soft parton emissions with 
small transverse momenta, which affect the precise At data because of constraints imposed on the 
transverse momentum pxi of the decay charged lepton. The At data require Pti to be above 20-25 
GeV in order to suppress charged leptons from background processes that do not involve a W boson 
decay. In addition, the Run-2 At data are organized into bins of Pti-, e.</., 25-35 GeV and 35-45 
GeV in order to better probe the x dependence of d(x, Q)/u(x, Q) [9]. While such binning amplifies 
the sensitivity of the At data to the PDFs, it also makes it dependent on the shape of the prt 
distribution near the Jacobian peak at pxt ~ M\y/2 « 40 GeV or, equivalently, to the transverse 
momentum (Qt) distribution of W bosons at Qt < 20 GeV, where large logarithms \x\(Qt/Q) 
dominate the cross section. A calculation that evaluates these logarithms to all orders in a s [25, 26, 
27], in addition to including the leading NNLO corrections [29, 28], results in somewhat different 



Figure 3: Comparison of the CT10W and CTEQ6.6 predictions with the D0 Run-II data for the 
electron charge asymmetry A e (y e ) for an integrated luminosity of 0.75 fb" 1 [10] 




Figure 4: The d/u ratio for CT10 (left) and CT10W (right) versus that for CTEQ6.6, at scale 
At = 85 GeV. 
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predictions for Ai than (N)NLO calculations without resummation, like those implemented in the 
other available codes [30, 32, 33, 34]. 

In CT10 fits, the QCD radiative contributions to Ae(yg) are implemented to the next-to- next- 
to-leading accuracy in Qt logarithms and NLO accuracy in the QCD coupling strength using 
the program ResBos that realizes the approach of Refs. [25, 26, 27]. The resummed differential 
distributions for da/dQx and da/dpxe both agree well with the data, in contrast to the fixed-order 
results. We thus expect that the resummed predictions for Ag implemented in the CT10 fit are 
more reliable as well. 

The magnitude of differences between the NLO and resummed predictions is illustrated by 
Fig. 2, comparing the CDF Run-2 A e (y e ) data [9] with LO, NLO, and resummed NNLL-NLO 
predictions from ResBos. The NLO and resummed curves are clearly distinct in the bin 25 < pre < 
35 GeV, shown in the left panel, and some differences are also seen in the bin 35 < pre < 45 GeV, 
shown in the right panel. The shape of the NLO prediction in this comparison is not unique and 
depends on the phase space slicing parameter Q^ P that defines the size of the lowest Qt bin where 
the real and virtual NLO singularities are canceled [26] . In the current comparison, Q^ P = 3Ge V, 
but other values of Q T ep are equally possible and would lead to NLO predictions lying closer to, 
or further from, the shown resummed curve. Such variations due to Q T ep or factorization scale 
indicate that resummation effects are important and should be included in precise fits to Ag. 1 

Numerical predictions of the CT10W analysis. When included in the CTEQ global 
analysis, the CDF measurements of Ag in Tevatron Run-1 [37] and Run-2 [9] agree well with the 
other data sets constraining the d/u ratio, provided by deep inelastic scattering on proton and 



lr The NNLO contributions produce marginal modifications compared to the NNLL-NLO result included in the 
CT10 analysis. We examined these contributions by redoing the calculation for A^(y^) after adding the exact o? s 
correction for W bosons produced with non-zero transverse momentum, which captures a large part of the full 
NNLO effect. The changes were found to be small and inconsequential in the current fit. 



deuteron targets by the NMC [38] and BCDMS [39, 40] collaborations. However, the most precise 
Run-2 lepton Ag data by the D0 Collaboration [10, 11] run into disagreement with the NMC and 
BCDMS deuteron DIS data, and in addition, exhibit some tension among themselves. Because of 
these disagreements, two separate fits are produced: CT10, which does not contain the D0 W 
electron asymmetry sets, and CT10W, in which they are included using weight factors larger than 
1 to ensure an acceptable fit. We obtain x 2 /Npt = 91/45 = 2 in the CT10W fit to the Run-2 
Ag, which is a significant improvement compared to 200/45 = 4.4 obtained in CT10. Comparison 
of NNLL-NLO predictions based on CT10W PDFs with the electron charge asymmetry data is 
presented in Fig. 3. 

In the CT10W analysis, the inclusion of the Run-2 Ag data increases the slope of d{x)/u{x) at 
x between 0.1 and 0.5 and reduces its uncertainty, as compared to CTEQ6.6 and CT10. This is 
illustrated by Fig. 4, which shows uncertainty bands for the d/u ratio in CTEQ6.6, CT10, CT10W 
PDFs vs. the momentum fraction x at scale Q = 85 GeV . 

Fig. 5 shows the ratios rwz — &{pp — > W ± X)/a(pp — > Z°X) and r w + w - = cr(pp — > 
W + X) / a{pp — > W~X) of the rapidity distributions in W ± and Z boson production at the LHC, 
obtained using CTEQ6.6, CT10 and CT10W PDFs and divided by the predictions based on the 
CTEQ6.6M set. Here, the reduction of the uncertainty bands in the ratio of W + to W~ cross 
sections predicted based on the CT10W PDFs, as compared to CT10, is again evident. 

Comparison with other fits. The MSTW'08 [5] and NNPDF2.0 [36] groups have also 
explored the impact of the D0 Run-2 Ag data. While their conclusions are broadly compatible 
with ours, the details vary. For example, in the CT10W fit, we find that all three pxg bins of the 
electron Ag data and the second pxe bin of the fit can be combined with the other data sets, despite 
the remaining disagreement with the deuteron DIS experiments. NNPDF, on the other hand, finds 
that all bins of the muon Ag are compatible with the deuteron DIS and all other experiments, but 
incompatible with the second and third prg bins of the electron Ag. NNPDF compute their NLO 
Ag predictions using the DYNNLO code [34], which can deviate from the resummed NNLL+NLO 
predictions by ResBos used by CTEQ, as discussed above. 

The resolution of the Run-2 Ag puzzle thus seems to require consistent implementation of 
perturbative QCD calculations both in experimental analyses and PDF fits, including small-Q^ 
resummed contributions, when constraints on the lepton pxg are imposed. For greater precision 
it may be preferable to perform the small-Qr resummation for every component of the angular 
distribution of the decay lepton [41], in addition to the resummation for two dominant angular 
components that is currently implemented in ResBos. Resummation of the full dependence on the 
polar and azimuthal angles 9 and <fi of the lepton in the vector boson rest frame may be important 
especially in situations in which the experimental coverage is not uniform in all directions, or when 
there are gaps in the coverage. Normally the angular coverage varies with rapidity, notably at large 
rapidity, and this variation may well affect the lepton asymmetry, with different consequences in 
different experiments. 



Figure 5: CT10, CTIOW, and CTEQ6.6 PDF uncertainty bands for the ratios 
(da(W ± ) / dy) / {da(Z) / dy) (upper two subfigures) and (da(W + )/dy)/(da(W~)/dy) (lower two sub- 
figures), at the LHC energies 7 and 14 TeV. 
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CT10 predictions for collider observables 

Fig. 6 compares the NLO total cross sections, obtained using CT10 and CT10W PDFs, to those 
obtained using CTEQ6.6 PDFs, for some selected processes at the Tevatron Run-2 and the LHC 
at yfs = 7 TeV. For most of the cross sections, CT10 and CT10W sets provide similar predictions 
and uncertainties, which are also in good agreement with those from CTEQ6.6 (i.e., well within 
the PDF uncertainty band). At the LHC, the PDF uncertainties in CT10 and CT10W predictions 
for some processes are larger than those in the counterpart CTEQ6.6 predictions, reflecting the 
changes in the framework of the fit discussed in the next paragraph. At the Tevatron, the CT10(W) 
PDF uncertainties tend to be about the same as those for CTEQ6.6, with a notable exception of 
ti production cross sections, which have a smaller PDF uncertainty with the CT10W set, because 
of stricter constraints on the up- and down-quark PDFs at the relevant x values. 

New PDF parametrizations; advancements in statistical analysis 

The CT10 global analysis implements several new features which were not available in the previ- 
ous studies. The systematic uncertainty associated with the overall normalization factor in each 
of the data sets is handled in the same manner that all other systematic error parameters are 
handled. The best-fit values of the normalizations are found algebraically, and their variations 
are included in the final estimate of the PDF uncertainties. More flexible parametrizations are 
assumed for the gluon, d-quark, and strange quark PDFs at the initial scale 1.3 GeV, to reduce 
biases in predictions in kinematical regions where the constraints from the data are weak. Finally, 
a new statistical procedure is introduced to guarantee the agreement of the fits at 90% C.L. with 
all included experiments, for any PDF eigenvector set produced by the error analysis. This is 
realized by adding an extra contribution to the total y 2 , which guarantees the quality of fit to each 
individual data set and halts the displacement along any eigenvector early, if necessary, to prevent 
one or more individual data sets from being badly described. The old procedure to enforce the 
90% C.L. agreement with all experiments in the CTEQ6 family of fits, by artificially increasing 
statistical weights of \ 2 contributions from those experiments that may be fitted poorly by some 
PDF eigenvector sets, is phased out by this more efficient method. As a result of these changes, the 
CT10/CT10W PDF uncertainty may be smaller or larger than the CTEQ6.6 uncertainty, depend- 
ing on whether the improved constraints from the data outweigh the increased uncertainty due to 
the relaxed PDF parametrizations and variations of normalizations during the determination of 
PDF eigenvector sets. 

Dijet invariant mass distributions 

Fig.7 compares NLO predictions based on CTEQ6.6, CT10, and CT10W PDFs with the data on the 
dijet invariant mass distribution da/dMjj reported recently by the D0 Collaboration [42]. These 
data are not included in the CT10 fit, but they are sensitive to the same scattering subprocesses, 
and include the same events, as the single-inclusive jet data constraining the gluon PDF in the 
CT10 fit. The cross sections are normalized to the theory prediction based on the central CT10 
PDF set. The statistical and systematic errors of the D0 data are added in quadrature. The 
renormalization and factorization scales in the theoretical predictions are set equal to a half of 
the average px of the jets, (pr)/2, consistent with the scale used in the single-jet cross sections 
when determining the gluon PDFs. With this choice of the scale, all three PDF sets agree with 
most of the data points within the PDF uncertainty. There appears to be a systematic excess 
of theory over the data, but its magnitude strongly depends on the assumed factorization scale. 
For example, it can be much worse if a different scale is taken, such as fi = (pr) assumed in Fig. 



Figure 6: Ratios of NLO total cross sections obtained using CT10 and CT10W to those using 
CTEQ6.6M PDFs, in various scattering processes at the Tevatron Run-II and LHC. 
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Figure 7: Comparison of D0 Run-II data for dijet invariant mass distributions [42] with NLO 
theoretical predictions and their PDF uncertainties for CTEQ6.6 (black), CT10 (red) and CT10W 
(blue) PDFs. The cross sections are normalized to theoretical predictions based on the best-fit 
CT10 set, designated as CT10.00. 
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2 of the D0 paper [42]. We conclude that the CT10 PDFs are reasonably compatible with the 
D0 dijet data within the present theoretical uncertainties, although an overall systematic shift, of 
order of the systematic shifts observed in the CT09 study of the single-inclusive jet distributions 
[1], would further improve the agreement, once the full correlated systematic errors of the dijet 
data become available. 



Uncertainty due to a s in CTEQ6.6 and CT10 PDF analyses 

Many calculations for collider processes (e.g., production of ti pairs and Standard Model Higgs 
bosons) require to evaluate two leading theoretical uncertainties, due to the choice of the PDF 
parametrization at the initial scale, and the value of the strong coupling constant a s (Mz) assumed 
in the hard cross sections and the PDFs. These uncertainties can be comparable in size, and their 
interplay, or correlation, may be important. In Ref. [2], we examine the a s dependence of CTEQ6.6 
PDFs [7] and show how the PDF-a s correlations are adequately captured by a simple calculation, 
without resorting to more elaborate methods proposed in other studies [43]. 

At the beginning of the PDF fit, one decides which data determine the a s (Mz) value and its 
uncertainty. CTEQ best-fit PDFs and their parametrization uncertainties are usually determined 
for a constant value of a s (Mz) that is close to its latest world-average central value; e.g., a s (Mz) = 
0.118 ±0.002 assumed in Ref. [2]. This input value of a s (Mz) can be viewed as an additional data 
point that summarizes world constraints on a 8 , mostly determined by precise experiments that 
are not included in the global fit (notably, LEP event shapes and t and quarkonium decays). In 
Ref. [2], we explore a more general procedure, in which the world-average data point on a s (M z ) is 
included in the fit in addition to the usual hadronic scattering data. A theoretical parameter for 
a s (Mz) is varied in this fit; its output value and uncertainty are determined by all input data. We 
find that the output value of a s (Mz) = 0.118 ± 0.0019 obtained in this way essentially coincides 
with its input value a s (Mz) = 0.118 ± 0.002. If the input value is not included, the output 
uncertainty on a s is increased significantly, to a s (Mz) = 0.118 ± 0.005. This indicates that the 
hadronic scattering data included in the fit imposes significantly weaker constraints on a s than the 
other experiments contributing to the world-average value. 

The Hessian PDF eigenvector sets returned by such floating-a s fit can be used to estimate 
the correlation between the PDF parameters and a s (Mz). However, each of these eigenvector 
sets is inconveniently associated with its own value of a s (Mz). Instead, one can apply a simpler 
procedure, in which all eigenvector sets except two are determined for the best-fit value of a s (Mz). 
These eigenvector sets provide the usual PDF uncertainty for a fixed a s (Mz). Separately, the 
uncertainty in the PDFs induced by the uncertainty in a s (Mz) is assessed, by producing two 
alternative PDF fits for the a s (Mz) values at the lower and upper ends of the a s uncertainty 
interval (i.e., a s (Mz) — 0.116 and 0.120). These PDF and a s uncertainties are then added in 
quadrature to obtain the total uncertainty. 

This procedure is valid both formally and numerically. It is based on a theorem that is applicable 
within the quadratic approximation for the log-likelihood function \ 2 in the vicinity of the best fit. 
The proof of the theorem, as well as a numerical demonstration of the equivalence of the addition 
in quadrature to the full estimation of the PDF+a s uncertainty based on the Hessian method, are 
given in Ref. [2]. The series of best- fit PDFs for a s (Mz) values in the interval 0.113-0.123, needed 
to evaluate the combined PDF+a s uncertainty in any application, are made available both for 
CTEQ6.6 and CT10 PDFs. 



Figure 8: Left: Ratio of g(x,Q) in SUSY fits with a fixed a s (M z ) = 0.118 and CT10 fit. Right: A x 2 in 
2004 and 2010 SUSY fits vs. gluino mass. 




Constraints on new physics from a global PDF analysis 

Besides providing the PDFs and their uncertainties, the global QCD analysis can establish bounds 
on masses of hypothetical particles beyond the standard model (BSM), for example, relatively 
light color-octet Majorana fermions that contribute to strong interaction processes with the same 
coupling strength as the gluons. Gluinos of a supersymmetric (SUSY) origin serve as an example 
of such fermions; but other models can introduce them as well. Constraints on the "gluinos" are 
often imposed in the context of a specific BSM model for their production and decay, which helps 
to rule out "gluinos" with masses up to a few hundred GeV. But, if no model-specific assumptions 
are imposed, much lighter gluinos are allowed: as light as w 12 — 15 GeV according to the 2004 
global fit based on the CTEQ6 set data [44], or 6-51 GeV according to the NLO/resummation 
analyses of e + e~ hadroproduction at LEP [45, 46, 47]. Note that the limit based on the global fit 
is not affected by (potentially important) theoretical uncertainties in the LEP analyses, associated 
with nonperturbative and matching effects in the resummation techniques that they employ. 

In Ref. [48], we improve the earlier limits [44] on relatively light gluinos based on an extended 
CT10 fit with added gluino scattering contributions. In this study, an independent PDF describing 
gluinos, and one-loop splitting functions describing interactions of gluinos with quarks and gluons, 
are introduced in the DGLAP equation. Two-loop gluino contributions are included in the renor- 
malization group equation for the running of a s . Cross sections for single-inclusive jet production 
are modified to include hard matrix elements for 2 — > 2 processes involving gluinos, with full 
dependence on the gluino mass evaluated in the general-mass factorization scheme. These mod- 
ifications capture the essential dependence on gluinos in inclusive processes studied in the CT10 
fits. Generally, the gluino contributions to inclusive observables are small, so that they can be 
evaluated at the one-loop level to achieve the same accuracy as the SM contributions evaluated at 
two loops. Squarks and other BSM particles are assumed to be heavier than a few hundred GeV 
and not included. 

Constraints on the gluino mass m g depend strongly on the value of a s (Mz) [44]. To reproduce 
the existing limits on a s {Q), we introduce two data points at Q = 5 GeV and Mz, representing a 
combination of measurements at low Q and Q w Mz, respectively. The low-Q bound on a s has 
been obtained by combining measurements of a s (Mz) in r and quarkonium decays. These low-Q 



measurements are not affected by gluinos heavier than 10 GeV. Gluino contributions to the high-Q 
data point are of the same order as the experimental uncertainties. 

With the a s data and latest hadronic data included, the 2010 SUSY PDF fits reduce the allowed 
range of gluino masses, as compared to the 2004 fits [44]. For example, Fig. 8 illustrates SUSY 
fits with a constant QCD coupling strength, a s (Mz) = 0.118, for gluino masses rrig shown in the 
figure. The left subfigure compares the gluon PDF obtained in SUSY fits with rrig = 20 and 50 
GeV (solid lines) to the CT10 error band (corresponding to rrig — oo). It is clear that too light 
gluinos distort the shape of the CT10 gluon PDF to an unacceptable level. The right subfigure 
shows the differences A\ 2 = x 2 ( a s, Trig) — Xctio- One can see that A% 2 > 100 for m g < 25 GeV, 
suggesting that gluinos lighter than 25 GeV are excluded at about 90% C.L., for a s {Mz) — 0.118. 
This improves the 2004 constraint [44], m g > 12 GeV for a s (M z ) = 0.118, by a factor of two. 
Similarly, the 2004 limit on rrig for a free a s , which allowed rrig = 1 GeV for a s {Mz) = 0.135, is 
increased to rrig > 13 GeV for any a s in the 2010 fit. 

Gluinos with mass about 50 GeV remain allowed both by the global fits and LEP data analysis. 
Such light gluinos may alter cross sections for jet production and other LHC processes [48]. We 
provide tables of PDFs with contributions of gluinos in this mass range to explore phenomenological 
implications. 




PDFs for leading-order showering programs 

Monte Carlo event generators, especially, the most mature leading-order generators, play a criti- 
cal role in all stages of modern particle physics. Neither conventional LO PDFs, nor NLO PDFs 
produce satisfactory results when implemented in the showering programs. In Ref. [4], we modi- 
fied the usual leading-order global analysis to find optimized PDFs for leading-order Monte-Carlo 
(LO-MC) simulations at the LHC. Besides the usual constraints from the existing hard-scattering 
experimental data, the joint input of this analysis incorporates pseudodata points for cross sections 
of W, Z, ti, SM Higgs, and gg — > bb' production at the LHC, as predicted by NLO QCD theory. 
The PDFs resulting from this analysis are not strictly at the leading order: they include some 
information about key LHC processes evaluated at NLO. Event generators, including "LO event 
generators", have some elements of higher-order contributions and, in this sense, are not at the 
stated order in the QCD coupling. We can use their flexibility to find the LO-MC PDFs that 
better reproduce the benchmark NLO cross sections, when combined with LO matrix elements in 
a fixed-order calculation or a LO event generator. 

To examine the available possibilities, we provide three representative LO-MC PDF sets, desig- 
nated as CT09MCS, CT09MC1, and CT09MC2. These PDFs realize different strategies for bring- 
ing the LO predictions closer to NLO. Their differences stem from varying assumptions about the 
running of a s (evaluated at one or two loops), factorization scales in the LHC pseudodata cross 
sections (fixed or fitted), and the momentum sum rule imposed on the PDFs (exact or relaxed by 
10-15% [49]). An example of the LO-MC PDFs in action is shown in Fig. 9, which compares cross 
sections for W + boson and Standard Model Higgs boson rapidity distributions at the LHC, ob- 
tained with LO matrix elements and CT09MC2 and MRST20071omod PDFs [49], and at NLO with 
CTEQ6.6. In the W + production case (upper left subfigure), the CT09MC2 calculation closely 
reproduces both the normalization and shape of the NLO cross section at all three LHC energies, 
while the MRST20071omod prediction differs from NLO in normalization at ^/s = 7 TeV, and 
both in normalization and shape at 10 and 14 TeV. For Higgs production (upper right subfigure), 
both CT09MC2 and MRST20071omod predictions provide almost identical distributions, which 
are smaller than the NLO prediction by a nearly constant normalization factor. This difference 
with NLO reflects especially large virtual corrections present in Higgs production cross sections, 
which cannot be completely compensated by an increase in the LO gluon density. However, since 
the average normalization factors K for each pseudodata process are also known from the fit (and 
published in our paper), end users can multiply the LO-MC cross sections for Higgs production 
and other pseudodata processes by these K-factors to better approximate the NLO cross sections 
(cf. the lower left subfigure). If an LHC process is not included as the pseudodata, comparison of 
LO predictions based on several LO-MC sets may still provide a reasonable estimate of the NLO 
cross section, as illustrated by the cross section for SM Higgs boson production via vector boson 
fusion in the lower right subfigure. 

In summary, the CT10 and CT10W sets are based on the most up-to-date information about 
the PDFs available from global hadronic experiments. There are 26 free parameters in both new 
PDF sets; thus, there are 26 eigenvector directions and a total of 52 error PDFs for both CT10 
and CT10W. The CT10 and CT10W PDF error sets, along with the accompanying a s error sets, 
allow for a complete calculation of the combined PDF+a s uncertainties for any observable. 
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